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Abstract. High-resolution X-ray spectra of high- mass stars and low- 
mass T-Tauri stars obtained during the first year of the Chandra mission 
are providing important clues about the mechanisms which produce X- 
rays on very young stars. For C, Pup (04 If) and (" Ori (09.5 I), the broad, 
blue-shifted line profiles, line ratios, and derived temperature distribution 
suggest that the X-rays are produced throughout the wind via instability- 
driven wind shocks. For some less luminous OB stars, like 9^ Ori C 
(07 V) and r Sco (BO V), the line profiles are symmetric and narrower. 
The presence of time- variable emission and very high-temperature lines in 
0^ Ori C and r Sco suggest that magnetically confined wind shocks may 
be at work. The grating spectrum of the classical T-Tauri star TW Hya 
is remarkable because the forbidden-line emission of He-like Ne IX and 
O VII is very weak, implying that the X-ray emitting region is very dense, 
ne ~ 6 X 10^^ cm~^, or that the X-rays are produced very close to the 
ultraviolet hotspot at the base of an accretion funnel. ACIS light curves 
and spectra of flares and low-mass and high- mass young stellar objects in 
Orion and p Ophiuchus further suggest that extreme magnetic activity is 
a general property of many very young stars. 



1. Introduction 

One of the early surprises from the Einstein mission was the detection of strong 
X-ray emission from O stars (Harnden et al. 1979). While the wind-shock 
paradigm for X-ray emission from hot stars (Lucy and White 1980; Owocki, 
Castor, & Rybicki 1988; MacFarlane & Cassinelh 1989; Hilher et al. 1993; 
Feldmeier et al. 1995; Cohen et al. 1996; Owocki & Cohen 1999) has gained 
acceptance over the past two decades, there is increasing evidence that not 
all early-type stars have X-ray properties that can be explained by instability- 
generated wind shocks alone. This evidence includes very hot plasma on the 
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Figure 1. Hard-band (2.5-8.0 keV: black circles) and soft-band (0.5- 
2.5 keV gray squares) light curves of a flare on the T Tauri star (GY 
195) in the p Oph A cloud. The ACIS flare spectrum can be fit 
with a single-temperature 4.7 ±0.4 keV, low-Fe abundance (Z = 0.25) 
VMEKAL plasma with Na = 3.9 ± 0.1 x 10^2 cm'^. The hard-band 
emission rises quickly at the onset of the flare and again during sub- 
sequent reheating events seen in the hardness ratio time-series (lower 
panel) . 



young BO V star r Sco (Cohen, Cassinelli, & Waldron 1997), overall levels of 
X-ray emission that approach or even exceed reasonable estimates of the wind 
emission measure in early- and mid-B stars (Cohen, Cassinelli, & MacFarlane 
1997), and strong 15.4-day periodic X-ray emission on the central O star of the 
Orion Nebula, 9^ Ori C (Gagne et al. 1997). 

The lack of outer convection zones in hot stars, and the lack of any correla- 
tion between X-ray activity and rotation (Pallavicini et al. 1981) has been taken 
to mean that the X-rays from early-type stars are not produced by magnetic ac- 
tivity associated with a solar-type dynamo. However, it has been suggested that 
large-scale magnetic fields, in conjunction with radiation-driven winds, could 
play a role in the X-ray production on some of these early- type stars (e.g., Ba- 
bel & Montmerle 1997a). The strong and hard X-ray emission from 6^ Ori C, 
r Sco, as well as several other young hot stars (Corocoran et al. 1994; Yamauchi 
et al. et al. 1996; Schulz et al. 2000), suggest that some type of magnetic con- 
finement, possibly associated with youth, may lead to the efficient production 
of hard X-rays on some hot stars. If the magnetic fields on hot stars survive 
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the star-formation process, the transfer of angular momentum via a magnetized 
wind would lead to significant spin down. 

The discovery of strong variable X-ray emission from classical T-Tauri stars 
(CTTS) in Taurus- Auriga with Einstein (Gahm et al. 1980) was less surprising 
because TTSs were known to be relatively fast rotators, strong Ha sources, 
and to undergo strong white-light flares. Since then, large-amplitude, long- 
duration flares have been observed many times with Einstein, ROSAT, ASCA, 
and Chandra. The temperature of the flaring plasma rises from a quiescent value 
of 6-20 MK to a peak value of 50-100 MK in minutes to hours and decays on 
time scales of many hours. The central question is whether TTSs have very 
active solar-type coronae driven by a magnetic dynamo or whether flares are 
produced in some kind of a magnetosphere coupled to a disk. 

The X-ray flare spectra and light curves have been modeled in two ways: as 
a single heating event in a quasi-statically cooling magnetic loop and as a scries of 
heating events in a coronal active region. E.g., a giant ASCA flare on the weak- 
line T-Tauri star (WTTS) V 773 Tau showed a peak temperature of 100 MK 
and a decay time of 2.3 hours (Tsuboi et al. 1998). A quasi-static cooling loop 
approximation implies moderate electron densities like those observed in solar 
flares and very large loop lengths, / ^ R^. Favata, Micela, & Reale (2001) 
have modeled the same flare data hydro dynamically to show that such flares 
may be produced in much smaller loops. In this model, long-duration flares are 
a manifestation of continuous reheating events rather than a single flare in a 
large, low-density loop. 

ASCA flares on three young stellar objects (YSOs) in p Ophiuchus, YLW 15, 
WL 6, and Elias 29, showed evidence of rotational modulation, quasi-periodic 
flaring, and significant spectral time variability (Tsuboi et al. 2000). Montmerle 
et al. (2000) present a magnetosphere/disk interaction model to explain the 
quasi-periodic flares on the Class I YSO YLW 15. Again, Favata, et al. (2001) 
show that such flares may reflect reheating events in high-latitude active regions. 

The recent discovery of X-ray emission from deeply embedded submillime- 
ter/radio dust condensations in the Orion molecular cloud (Tsuboi et al. 2001) 
indicates that magnetic heating begins at the earliest stages of star formation. 
It remains to be seen whether the large-scale magnetic fields which are believed 
to regulate accretion, mass-loss, and bipolar outflows are related to the recon- 
nection events which lead to flares and X-ray activity. 

Disk-photosphere coupling via a magnetic field was first suggested by Konigl 
(1991) (inspired by the model of Ghosh k. Lamb 1979) to explain the slow 
rotation of CTTSs. A recent study of low-mass 1-30 Myr stars in the Orion 
Nebula Cluster and NCC 2264 (Rcbull et al. 2001) found that PMS stars evolve 
with nearly constant angular velocity as they contract down their convective 
tracks. This means that the mechanism which regulates the spin of a CTTS 
may still be at work on WTTSs long after a dusty disk has been dissipated. If 
magnetic coupling still operates long into the WTTS phase, then a gas (plasma) 
disk must exist within a few stellar radii of the photosphere. 

In this paper we will discuss recent Chandra results which address magnetic 
activity on young high-mass and low-mass stars. §2 summarizes results based on 
ASCA and Chandra light curves and CCD spectra of YSO flares in p Ophiuchus. 
In §3 we describe the X-ray grating spectrum of the CTTS TW Hya and use the 
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forbidden lines of Ne IX and O VII to suggest a location for the X-ray emitting 
plasma. In §4 we describe grating spectra of four hot stars that, in some cases, 
provide strong evidence for magnetic confinement. In §5 we present 2-D MHD 
simulations of magnetically confined wind shocks. Finally, we conclude that 
magnetic fields, disks, and mass-loss may play a central role in the production 
of X-rays on very young high-mass and low-mass stars. 

2. Young Stellar Objects in p Ophiuchus 

The p Ophiuchus A cloud was observed continuously for 96 ks with ACIS-I 
on 2000 May 15. An RGB image of the fiterered event data with rcd=soft 
(0.5-1.5 keV), green=medium (1.5-2.5 keV), and blue=hard (2.5-8 kcV) shows 
76 X-ray sources. Highly absorbed YSOs and background AGN seen through 
25-75 mag of visual absorption appear as blue sources in this image. Twenty- 
two background AGN have been identified as faint sources with high 2-keV 
hardness ratios and low Kolmogorov-Smirnov variability statistics. By dividing 
the event data into 31 time bins and removing flickering background pixels, we 
have produced a color movie illustrating what Montmerle et al. (1983) described 
as an "X-ray Christmas Tree" in p Oph. 

In Figure 1, we present ACIS medium- and hard-band light curves of the 
low-mass YSO GY 195 in p Oph A. The lower panel is the time-series of 3-keV 
hardness ratios for GY 195. This flare is interesting because of the deep dip in 
the hard-band (3 — 8 keV) light curve approximately 48 ks after the start of the 
observation. Initially, we interpreted this as evidence of self-eclipsing. However, 
the soft-band (0.5 — 3 keV) light curve shows a much smaller dip than expected. 
Since the emission-measure distribution at the peak of a flare is dominated by 
very hot plasma (T > 50 MK), those times of the when the hard emission rises 
sharply probably correspond to reheating events. We see a few small flare events 
before the major peak at 41 ks, a rapid decay, and a series of reheating events at 
approximately 52, 62, 83, and 93 ks. We have yet to perform a detailed spectral 
time-series analysis on these data, but our initial XSPEC analysis suggests a 
very hot, low-abundance coronal plasma, and high column density consistent 
with other embedded T Tauri stars in p Oph. 

The p Ophiuchus cloud contains two optically visible stars: the young, mag- 
netic B3 star SI Oph and the WTTS DoAr 21. During the Chandra observation, 
DoAr 21 was in the decay phase of a very large flare. Although no major flares 
were seen on SI, the light curve exhibits low-amplitude, long-term and short- 
term variability. The ACIS spectrum of this high-mass YSO is reminiscent of 
low-mass YSOs: iVn = 4.1±0.5 x 10^2 cm-^ and kT = 3.1±0.5 keV. Simultane- 
ous VLA observations at 2 and 6 cm show variable, non-thermal radio emission, 
confirming its magnetic origin, first reported by Andre et al. (1993). Si's radio 
and X-ray properties make it a good candidate for magnetically confined wind 
shocks (see §5). 

3. The Classical T Tauri Star TW Hya 

With the launch of Chandra and Newton-XMM, high-resolution grating spec- 
troscopy has become an important new tool for studying X-ray activity on young 
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Figure 2. The MEG speetrum (solid lines) of TW Hya around the rif 
Ne IX and O VII triplets (Kastner et al. 2001). The observed spectrum 
is overlaid with a variable-abundance, differential emission-measure 
model (dashed curve) that best-fits temperature-sensitive lines. Note 
the discrepant intercombination and forbidden lines. 



stars. The He-like resonance (r), intercombination (i), and forbidden (f) lines of 
coronal sources like AB Dor (dotted line in Fig. 2) are useful density diagnostics 
because the collisional excitation rate from the metastable ^Si state depends on 
density (see Figure 1 of Gabriel & Jordan 1969). The critical density is reached 
when the collisional excitation rate equals the forbidden '^Si ^Sq photodecay 
rate. For AB Dor, the forbidden-line emission is strong; the S, Si, Mg, and 
Ne f /i ratios are close to their low-density limits. The O VII f /i ratio yields 
logne = 10.68 ±0.33 cm-^. 

For TW HYa, a classical T Tauri star, the / lines of Ne IX (seen in Fig. 1) 
and O VII arc remarkably suppressed. Kastner et al. (2001) measure f/i = 
0.44 zb 0.13 for Ne IX indicating logne = 12.75 ± 0.15 cm"^. This is the highest 
coronal density measured so far using f/i ratios, exceeding that of AB Dor by 
two orders of magnitude, possibly signaling a fundamental difference between 
coronal stars and T Tauri stars. Kastner et al. go further, suggesting that the 
bulk of the X-ray emission on TW Hya is generated at the base of the accretion 
funnel which couples the accretion disk to the stellar photosphere. 

We point out, however, that ultraviolet emission can depopulate the metastable 
^Si state via photoexcitation to the ^P2io state (see Kahn et al. 2001). The 
forbidden line is suppressed when the '^Si '^P photoexcitation rate RpE ap- 
proaches the ^Si —>■ ^So photodecay rate Rrb- For O VII and Ne IX photoex- 
citation occurs in the ultraviolet at 1616 and 1246A. The photoexcitation rate 
for the metastable state of Ne IX is Rpe = F^^f, where / = 0.0949 (Cann 
& Thakkar 1992). Rpe = Rrd = 1-09 x 10"^ s^^ (Drake 1971) imphes a critical 
flux = 4.33 X 10^ s"^ cm-2 Hz'^ TW Hya's flux at 1245.8 A from the lUE 
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Figure 3. MEG spectra (solid gray lines) of C Pup (upper left panel) 
and 6^ Ori C (lower panel). The MEG spectrum of AB Dor (KO V) is 
shown (dashed line) for comparison. As expected for a luminous hot 
star with instability driven wind shocks, the Ne X resonance line of 
C Pup is blueshifted, broad, and asymmetric. The Ne X line of 6^ Ori C 
(upper right panel) is symmetric about line center and its non-thermal, 
non-instrumental width is f ~ 410 km s^^. Note the strong Fe XXI 
line. C Ori (lower left panel) shows a broad Ne X line. The r Sco lines 
(lower right panel) are narrow, like those of the coronal source AB Dor 
(dashed line). 



archive is f\ « 3.5 x 10~^^ ergs s~^ cm~^ A . We consider the simplest geo- 
metrical model: a star at a distance D (56.4 pc for TW Hya) with X-ray and 
UV sources separated by a distance d. In this case, Fj^ = {D/d)'^/^. Solving for 
the UV/X-ray separation at which i?PE = -RpD, we find d « 2.5 x 10^ cm, i.e., 
d pa QMRq. 

On TW Hya, the O VII forbidden line is strongly suppressed. Assuming 
RpB = 1.04 X 10^ s~^ and / = 0.0742 implies a critical incident flux Fi, = 
5.28 X 10^ s-i cm-2 Hz'^ TW Hya's lUE flux at 1616.5 A is fx ^ 3.0 x 
10~^^ ergs s~^ cm~^ A ^, implying a critical UV/X-ray separation d 1.1 x 
10^° cm. Since the O VII / line is entirely suppressed, d < 0.16i?Q. Thus, the 
Ne IX and O VII f /i ratios may be suppressed via ultraviolet photoexcitation, 
provided the X-rays are produced close to the UV continuum source. 

We note that T Tauri stars have strong UV emission lines (H I, N V, C IV, 
Si IV) produced by optically thin 5 x 10^ — 2 x 10^ K plasma. The continuum 
emission (which depopulates the state) is not produced by an optically 
thin transition region. If this thick material is distributed over much of the 
photosphere, the UV continuum would be too diluted to suppress the O VII 
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and Ne IX forbidden lines. The forbidden lines are suppressed if the X-rays are 
emitted at the base of the accretion funnel directly above the optically thick UV 
hotspot. 

An alternate hypothesis is that X-ray emitting magnetic loops are dis- 
tributed throughout an active corona. In this case, the f /i ratios are unaffected 
by UV flux and logrie « 12.75 cm~^, the total particle density n = 1.77nc, the 
mean temperature in the loop is at the peak of the emission-measure distribu- 
tion T ~ 3.0 X 10^ K, and the total gas pressure in the loop is P = nkT = 
4.12 X 10^ dynes cm~^. The mean magnetic field strength required to confine 
such a loop is > SirnkT ^ B > 320 G. Direct measurements of magnetic 
field strengths and filling factors on other T Tauri star indicate very strong 1- 
4 kG fields concentrated in cool spots covering a large fraction of the photosphere 
(Johns-Krull, Valenti, Sz Koresko 1999; Guenther et al. 1999). Even allowing 
for the geometric reduction in the magnetic intensity as loops expand into the 
corona, coronal loops on T Tauri stars can sufficiently confine high-density coro- 
nal plasma. 

In conclusion, X-ray and far-UV spectra of T Tauri stars may provide a key 
test of accretion-powered X-ray activity. If the forbidden lines are suppressed 
by photoexcitation, Si XIII, Mg IX, Ne IX, and O VII f/i ratios will yield 
self-consistent estimates of the UV/X-ray source separation. 

4. X-ray Grating Spectra of OB stars 

The far-ultraviolet flux from OB stars is 10^ — 10^ times brighter than on T 

Tauri stars like TW Hya. As a result, O VII and Ne IX f-linc emission is 
essentially absent. The ^Si transitions of Mg XI and Si XIII at 1024.4 A 

and 901.8 A, respectively, are used to estimate the average distance bewteen the 
photosphere and the X-ray emitting region. Kahn et al. (2001) and Cassinelli 
et al. (2001) have shown for ( Pup (04 If) that the visible X-rays are emitted 
1 — from the photosphere where the overlying wind becomes optically thin 
at the wavelengths of Mg XI, Si XIII, and S XV. These results indicate that 
X-ray emitting regions are distributed throughout the wind as predicted by 
instability-driven wind shock theory. ( Pup's X-ray line profiles provide further 
confirmation of this picture. In the upper left panel of Figure 3, the 12.1 A NeX 
line of Pup is blueshifted, broad, and asymmetric because the optically thick 
wind absorbs most of the X-rays from receding shocked material on the far side 
of the star. The challenge for wind shock models is to produce strong enough 
shocks at the base of the wind to explain the Si XIII and S XV emission. 

While high-resolution spectra of ^ Pup support the distributed wind-shock 
hypothesis, HETG spectra of three other hot stars are more difficult to explain. 
In Fig. 3, the Ne X line of 9^ Ori C (07 V) and r Sco (BO V) do not show 
obvious line asymmetries. In their analysis of the 0^ Ori C spectrum, Schulz et 
al. (2000) find a range of line widths (300 — 800 km s~^) and no systematic line 
shifts. Schulz et al. (2000) state that these results are consistent with instability- 
driven wind shocks. However, Owocki &; Cohen (2001) have shown that the 
substantial amount of emission longward of line center is difficult to reconcile 
with the expected attenuation by the wind. Based on f/i ratios, Waldron & 
Cassinelli (2001) find that the hottest plasma around Ori must lie less than 
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IRi, from the photosphere, concluding that wind shocks and some magnetic 
confinement of turbulent hot plasma may be required. The r Sco lines are 
narrow and may not be produced in outflowing material many stellar radii from 

the photosphere. 

The 10-12 A spectra of C, Pup, 6^ Ori C, and r Sco also exhibit striking 
differences. C Pup (see Figure 3 of Cassinelli et al. 2001 and Figure 5 of Kahn 
et al. 2001) shows broad, asymmetric, blue-shifted lines of Ne IX and Fe XVII, 
indicating little emission measure above 10 MK. The 10-12 A MEG spectra of 
9^ Ori C and r Sco in Figure 4 show very strong, narrow, symmetric Fe XXIII 
and Fe XXIV lines suggesting very hot (T > 16 MK), solar-abundance plasma 
(Schulz ct al. 2000). To date, instability- driven wind shock models cannot 
produce such high-temperature shocks, especially not on OB stars like 0^ Ori C 
and r Sco with relatively low mass-loss rates and terminal wind speeds. 



5. Magnetically Confined Wind Shocks 

The leading theoretical picture for X-ray production in magnetized hot stars 
is the magnetically confined wind shock (MOWS) model proposed by Babel & 
Montmerle (1997a) and applied by the same authors to 9^ Ori C (1997b). In this 
model, an ionized, radiation-driven wind is confined to fiow along a magnetic 
dipole field and is thus channeled, at low latitudes, into the magnetic equatorial 
plane. In the magnetic equator, the streams from the two hemispheres collide 
and are confined by the dipole field, producing a standing shock and X-ray 
emission above and below a disk-like feature above the magnetic equator. 

This model succeeds in producing relatively hard X-rays (due to the nearly 
head-on collision of two fast streams), high levels of emission (due to the con- 
finement at relatively high densities of the shock- heated plasma) , and rotational 
modulation of the X-rays, like that observed in 0^ Ori C, if the dipole field is 
tilted with respect to the rotational axis. This oblique magnetic rotator model 
has been used to explain many of 9^ Ori C's curious behaviors: periodic Ha 
and He II emission and periodic C IV and Si IV P-Cygni profiles (Stahl et al. 
1993; Walborn & Nichols 1994; Stahl et al. 1996; Reiners et al. 2000). Based on 
circular spectro-polarimetric monitoring of its Balmer lines, Jean-Francois Do- 
nati and collaborators have recently measured a polar magnetic field strength 
B 1.1 kG (2001, private communication) at phase 0.0. This implies that 
X-ray, Ha and He II emission is greatest when looking down on the magnetic 
pole (phase 0.0). Emission is weakest when looking at the magnetic equator 
(phase 0.5), presumably because hot plasma just above the magnetic equator is 
eclipsed by the photosphere. We note that Stahl et al. and Babel & Montmerle 
proposed the opposite magnetic geometry. 

The model of Babel & Montmerle, which assumes a rigid magnetic geometry, 
has been extended to allow the magnetic field morphology to be modified by the 
wind flow in a self-consistent manner (Ud-Doula &: Owocki 2001, in preparation). 
Assuming an initially dipolar field and spherical mass-loss, the 2-D MHD code 
simulates the time evolution of density, temperature, and velocity throughout 
the wind. A given set of field strength and wind parameters defines a confinement 
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Figure 4. Upper panel: the MEG 10-12 A spectrum of 6^ Ori C 
shows strong, moderately broad, symmetric Fe XX-Fe XXIV lines in- 
dicating substantial emission measure at or above 16 MK. Lower panel: 
the T Sco (BO V) spectrum shows narrow Fe XX-Fe XXIV lines similar 
to those seen on the active corona of AB Dor (Kl V) (dashed lines). 
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Figure 5. 2-D MHD simulations of an initially dipolar magnetic field 
distorted by radial mass loss on ^ Pup. Density (upper panels) and non- 
radial velocity contours (lower panels) at the end of two realizations: 
moderate magnetic confinement {rj = VTO, t = 450 ks) on the left and 
high confinement {rj = 10, t = 295 ks) on the right. With sufficient 
confinement, material is accelerated to approximately ±1000 km 
creating shocks in dense, unstable regions near the magnetic equator 
only 0.1 — IRi, above the photosphere. The recent measurement of a 
1.1-kG polar field on 9^ Ori C implies 77 13 as on the right. The 
resulting dense, high-speed shocks may explain the broad, symmetric, 
high-temperature lines seen on 9^ Ori C. 
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parameter, 

where Bq is the surface equatorial field strength (approximately half the polar 
field strength) . 77 3> 1 implies the wind is confined by the field and 77 <C 1 implies 
the field is completely stretched out by the wind. 

In Figure 5, we show a pair of C Pup simulations for rj = v^To and r] = 10. 
For 9^ Ori C, Bq ~ 550 G and rj ~ 13, implying high-speed shocks, T ^ 10^ K, 

and rie <; lO^'' cm^^. As of this writing, we do not know if the wind continuum 
optical depth is low enough to view the entire volume of X-ray emitting plasma. 
Nonetheless, the 2-D simulations suggest that a moderately strong dipolar field 
can provide significant confinement and shock heating. To realistically test the 
MCWS idea, synthetic X-ray spectra based on 3-D MHD simulations of an 
oblique magnetic rotator are needed for a range of mass-loss rates, magnetic 
field geometries, and viewing angles. 

6. Discussion 

The velocity profiles, f /i ratios, and high plasma temperatures derived from 
Chandra spectra of OB stars suggest that magnetic fields may play an impor- 
tant role in the production of X-rays on young, massive stars. Time-resolved, 
high-resolution X-ray spectra of very young stars may help uncover the geome- 
try of the underlying field. It is interesting to note that SI Oph, 9^ Ori C, and 
T Sco, which show the most striking signs of magnetic activity, are all young 
and associated with star- forming regions: p Oph 1 Myr, Andre et al. 1988), 
Orion Nebula 1 Myr, Hillenbrand 1997), and Sco-Cen 10 Myr, Kihan 
1994). Although a different kind of magnetic activity produces very hot, ener- 
getic X-ray flares seen on various classes of low-mass YSOs, both high-mass and 
low-mass stars undergo a brief episode of rapid accretion and mass-loss regu- 
lated by large-scale fields and disks. The evolution of the field and disk affects 
angular momentum evolution, planet formation, and the high-energy radiation 
environment of young planetary atmospheres. 
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